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Abstract: Here we employ both dynamic and geometric phase control 

of light to produce radially modulated vector-vortex modes, the natural 
modes of optical fibers. We then measure these modes using a vector modal 
decomposition set-up as well as a tomography measurement, the latter 
providing a degree of the non-separability of the vector states, akin to an 
entanglement measure for quantum states. We demonstrate the versatility of 
the approach by creating the natural modes of a step-index fiber, which are 
known to exhibit strong mode coupling, and measure the modal cross-talk 
and non-separability decay during propagation. Our approach will be 
useful in mode division multiplexing schemes for transport of classical and 
quantum states. 
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1. Introduction 

A key technique in achieving high-capacity optical communication is that of mode division 
multiplexing (MDM), which uses the spatial modes of light as information carriers m , and has 
been shown promise in both free-space 0|3) and fiber |4|[5). To date there has been work on 
MDM with orbital angular momentum (OAM) Bessel beams El and linearly polarized 
(LP) modes commonly used to describe fiber modes M- However, within free-space, atmo¬ 
spheric turbulence has been shown to negatively affect the strength of communication signals 
m, and so MDM in fibers remains topical despite the challenges. 

The propagation of scalar spatial modes through fibers (the LP modes) has been limited to 
only a few modes, due to intermodal coupling within the fiber 0- This coupling is caused by 
inner fiber imperfections, such as index inhomogeneities, core ellipticity and eccentricity, as 
well as external perturbations like bends or strain CD- Despite these difficulties, successful 
fiber communication demonstrations have made use of specially designed fibers Ifl2l[l3l . and 
more recently using cylindrical vector modes mna. These modes are close to the natural 
modes of optical fibers as they consist of spatially inhomogeneous states of polarization lfl6l 
H3, but typically lack the radial modulation of the pure fiber modes. 

In this paper, we exploit both the dynamic (optical path length) and geometric (Pancharat- 
namBerry) phase changes of light to generate the natural modes of fibers: radial, azimuthal 
and polarization control. We use a spatial light modulator to produce the radial modes before 
creating vector-vortex versions by geometric phase control M- We introduce a composite di¬ 
agnostic in the form of a vector mode decomposition that is mode specific 0, as well as a new 
tomography approach that returns a quantitative measure of the degree of non-separability or 
’’vectorness” of these modes ED. We use our set-up to create modes for a step-index fiber with 
known mode coupling 11201 and demonstrate that we are able to quantify the propagation of the 
natural modes of these fibers for the first time. We quantify the effect of fiber perturbations on 
vector modes, which we demonstrate experimentally. This approach will be useful in studies 
of MDM for increased bandwidth as well as in transport of high-dimensional quantum states 
down fibers. 

2. Concept and implementation 

The full vector wave equation for a step-index fiber is given by IfTD 

{vf + n 2 lr}et + Vt{efVt [ln(n 2 )]} =/3 2 e), (1) 

where n is the refractive index of the fiber, k = 2n/X is the wavevector, e t is the transverse 
electric field and p is the propagation constant of each vector mode solution. From Eq. 0, 
the first four higher-order vector solutions may be described as: transverse electric (TEoi), 
transverse magnetic (TMqi), hybrid electric odd (HE°j d ) and even (HEo™ 11 ). The two indices 



represent the number of half-wave patterns across the width and the height of the waveguide, 
respectively. The whole mode group is nearly degenerate, where the two HE 21 modes share the 
exact propagation constant, as shown in Fig.|T| As such, these modes are likely to couple during 
propagation within step-index fibers. Due to the cylindrical symmetry and step-like constant 


(a) fiber scheme 


(b) scalar inodes 


(c) vector modes 



Fig. 1. Schematic of the step-index fiber and its modes; (a) structure of the fiber composed 
of a protective coating, cladding and core where the light is guided; schematic of the mode 
patterns and spectrum for (b) the scalar approximation and (c) the natural vector modes 
with partialy resolved degeneracy of the first higher order mode group. 


refractive index distribution of these fibers, the radial field distribution is given by Bessel-type 
functions; 


E e P (r) 


J \e\(^f)/ J \e\( u ep) for r<a 

K\e\ {'^f)/ K \l\ i w e P ) for r>a 


( 2 ) 


Here, a is the core radius, Jg is the fth order Bessel function, K t is the /th order modified 
Bessel function with ui p and W( p are normalized propagation constants. In our experiment, we 
made use of a 30 mm long step-index fiber with a core radius a = 15 ffm and a numerical 
aperture NA = 0.08. For ease of alignment, we used a wavelength of 633 nm resulting in a total 
of 76 possible modes that can exist in the fiber. 

Owing to the radial distribution of the fiber, it is important to match the input field to the 
corresponding fiber mode group. As such, we exploit both the dynamic and geometric phase to 
generate the natural vector modes of the fiber, with our experimental setup shown in Fig. [2] The 



Fig. 2. Experimental setup used to generate and detect vector beams after propagating 
through a fiber by two separate measurement techniques: vector mode detection and to¬ 
mography. Insets depict beam profiles at different planes. SLM - spatial light modulator; 
HWP - half wave plate; QP - q-plate; MO - microscopic objective; L - lens; NPBS - non¬ 
polarizing beam splitter; LP - linear polarizer; QWP - quarter wave plate. 


first component of the experiment generated the natural fiber modes by radial, azimuthal and 



























polarization control, labeled mode generation in Fig. [2] Firstly, the fundamental fiber mode, 
described by Eq. ([2]) with £ = 0 and p = 1, was encoded onto a spatial light modulator (SLM) 
by complex amplitude modulation [22;]. Secondly, we made use of a q-plate l23li24l to prepare 
vector vortex modes. Polarization is coupled to OAM by the selection rules: \i,L) \t 4 Q,R) 
and | £,R) —> \£ — Q,L). The azimuthal charge introduced by the q-plate is Q = 2 q. In our experi¬ 
ment q= 1/2 and the incident beam generated on SLM 1 was horizontally polarized with £ — 0, 
which generated a TEoi or azimuthally polarized vortex mode consisting of a superposition of 
£ = ±1 modes. Thus, the field after the generation step can be described by 

|*P) = TE 0 i (r, 0) =E\\(r) (exp(/0) \R) — exp(—/0) |L)), (3) 

where the held singularity introduced by the q-plate, results in a variation in the radial compo¬ 
nent E( p as well as in the azimuthal (@ ) degree of freedom. Note the OAM and polarization are 
not separable in Eq. |T} and as such V F is reminiscent of a two-dimensional maximally entangled 
quantum state (qubit), a fact that we will exploit in our detection of the modes after the fiber. 

After propagation through the fiber the output modes were directed to the detection scheme 
comprising of two parts: a vector mode detection and a tomography measurement. The for¬ 
mer acts as a vector mode decomposition technique for the fiber mode group [O, which we 
implemented for the first four higher-order fiber modes. For simplicity, we describe the single 
channel as shown in Fig. [2] however in reality we implemeted a four-channel detection for each 
of the fiber modes, which were measured simultaneously. The output fiber mode was directed 
through a half-wave plate oriented at 0 rad followed by a q-plate (q = — 1 / 2 ); it can be shown 
that this combination forms an inner product measurement in an analogous manner to that used 
for scalar mode detection 12(71 . This pairing was followed by a linear polarizer, oriented to 
transmit horizontally polarized light. The final component of the decomposition measured the 
on-axis intensity on a camera after a Fourier lens L 3 [20]. When the input mode matches the 
filter for the detected mode, a bright on-axis intensity is observed, otherwise zero on-axis in¬ 
tensity is measured. This detection scheme therefore determines whether the fiber causes any 
intermodal coupling between the four first higher order fiber modes. 

The second analysis technique, the tomographic measurement borrows techniques from tra¬ 
ditionally quantum measurements EU and is used to determine the effect of the fiber on the 
degree of non-separability or vectorness of the output mode. As non-separablity is not unique to 
quantum mechanics, many of the tools used to measure quantum states are applicable to vector 
beams m. We now apply this tool for the analysis of the natural modes in fibers. Projective 
measures are first performed on the polarization state, while the OAM degree of freedom was 
measured using holograms encoded onto SLM 2. As SLMs are polarization sensitive, a lin¬ 
ear polarizer could not be used to measure the polarization states, as is commonly performed. 
Instead, a half-wave plate was inserted before SLM 2, and rotated to specific orientations to 
realize a filter for the linear polarization states: horizontal, vertical, diagonal and anti-diagonal. 
By inserting a quarter-wave plate, we were able to filter the two circularly polarized compo¬ 
nents, resulting in a total of six polarization measurements. Similarly, we created six different 
holograms on SLM 2 to represent the two pure OAM modes as well as four orientations of 
the superposition states: \£ = 1 ) + exp(/ 02 )|(? = — 1 ), for 02 = 0, zr/2, zr,3?r/2. A modal 
decomposition was performed for each polarization state by measuring the on-axis intensity on 
a camera situated after a Fourier lens L 4 . This tomographic method produces an over-complete 
set of 36 measurements, which can be used to minimize the Chi-square quantity and recon¬ 
struct the density matrix p = |yr) (i/r| of the state of the vector mode (e.g., that of Eq. (0). 
This approach is analogous to measurements of qubit states except here one degree of freedom 
is polarization and the other spatial. The degree of non-separability of any field can then be 
calculated from the density matrix. The concurrence of a state typically represents the degree 


of entanglement in a quantum entanglement experiment l26l . However, as the vector modes 
are non-separable in the polarization and OAM degrees of freedom, the concurrence provides 
a measure of the degree of non-separability of the vector mode. 

3. Experimental results and discussion 

Using the vector mode detection techniques described above, the results demonstrate inter- 
modal coupling within the fiber between the first higher fiber modes (as expected). Figure0a) 
shows an example of the vector decomposition results, where a TEoi mode was generated after 
the q-plate and the on-axis intensity was measured in the four channels. It is clear from the 
images that a signal was only detected in channel TEoi. The comparison of the decomposition 
performed before and after the fiber, as shown in Fig. [3jb) and (c), respectively, illustrates the 
intermodal coupling between the four modes within the fiber. 


Input mode TE 01 TM 0 i HE e 21 HE° 21 



(b) (c) 


Fig. 3. Example and results of the vector mode decomposition; (a) depicts the working 
principle by illumination with a TEqi mode. The bright spot in the TE channel on the 
optical axis, indicated by the green crosses, compared to zero intensity in the other channel 
denotes the high mode purity. Illustration of the modal transmission matrix through (b) 
free space, indicating the high quality of the vector beam generation, and (c) after fiber 
transmission, showing the appearance of mode coupling during propagation through the 
fiber. 

From the 36 tomographic measurements shown in Fig. [4} a), the qubit vector density matrix 
was calculated as shown in Fig. |4|b). We found the concurrence, a measure of the degree of 
non-separability, of our vectorial field to be 0.98 ± 0.01. A maximal vector field is represented 
by 1, while a value of 0 represents a purely scalar field. Thus, although our first measurement 
scheme reported intermodal coupling within the fiber, the output mode has maintained a high 
degree of non-separability. We also measured the degree of non-separability for a higher order 
radial mode group with f = 0 and p = 2 in Eq. ({2}. We found it to be lower than that of the first 
mode group with a value of 0.85 ± 0.01. Since this mode group shows a lower power fraction 
confinement within the core than the first mode group EH . this mode group was more strongly 
affected by fiber imperfections, thereby resulting in a stronger decay of non-separability. 

This technique enables the investigation of the effects of fiber perturbations on the degree of 













Fig. 4. Experimental results obtained from the full state tomography of the TEqi mode, 
for qualitative analysis by rotating polarizer (Medial), (a) Tomography matrix obtained 
by projection of the mode into six polarization and the corresponding six OAM states, 
measurement procedure (Media2). (b) The qubit vector density matrix showing the high 
degree of non-separability of the investigated mode. 


non-separability of the vector modes. Table Q] shows the effect of adding stress to a fiber. By 
stressing the fiber, the energy exchanged between the modes fluctuates in a sinusoidal manner. 
Depending on the level of stress, the coupling between the modes can be increased or decreased. 
As the coupling length is inversely proportional to the coupling strength, modes that couple 
weakly will have longer coupling lengths. Thus, increasing the stress on the fiber does not 
directly result in a decrease in the degree of non-separability and as such we measure a degree 
of non-separability of 0.43 and 0.73 for stresses of 0.5 N and 1 N, respectively. 


Table 1. Effects of stress on degree of non-separability 


Stress (N) 

Degree of non-separability 

0 

0.98 ±0.01 

0.5 

0.43 ± 0.01 

1.0 

0.73 ±0.01 


4. Conclusion 

Here we have demonstrated a complete generation and detection scheme for the natural modes 
in fibers. We use a spatial light modulator to create radially modulated modes through complex 
amplitude modulation, and then exploit geometric phase change to control the polarization and 
OAM of the modes. We have demonstrated a composite detection method with which to analyze 
the effects of fiber propagation on these modes, using step-index fibers as an example. This 
technique can be extended to investigate the fiber propagation of unknown vector beams Ez). 
As the non-separability of vector modes offer increased potential for optical communication in 
both the classical and quantum regimes and as such, the ability to quantify the degree of non¬ 
separability is a useful technique, and the first time it has been applied to fiber mode analysis. 
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